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The contribution of the strange-quark current to the electromagnetic form factors 
of the nucleon is studied using lattice QCD. The strange current matrix elements 
from our lattice calculation are analyzed in two different ways, the differential 
method used in an earlier work by Wilcox and a cumulative method which sums 
over all current insertion times. The preliminary results of our simulation indicate 
the importance of high statistics, and that consistent results between the varying 
analysis methods can be achieved. Although this simulation does not yet yield a 
number that can be compared to experiment, several criteria useful in assessing 
the robustness of a signal extracted from a noisy background are presented. 



1 Introduction 

An important theme of contemporary Hadron Physics is the role of nonva- 
lence degrees of freedom. In particular, the contribution of strange quarks to 
a variety of nucleon properties has been studied.EJ For nucleon form factors 
the contribution of the strange-quark current can be extracted using infor- 
mation obtained from parity violation in polarized- electron scatteringfl A 
number of experimental results have been reportedtra and new measurements 
are planned. As well, there are numerous calculations using a number of 
different approaches a 

In this work we focus on the calculation of the strange-quark current loop 
using lattice QCD. Calculations of these so-called disconnected current inser- 
tions have been reported previously.floEl However, these calculations differ in 
the method used to analyze the results of their lattice simulations and also 
differ in their conclusions. WilcoxB used a differential method to analyze the 
time dependence of the three-point function and no signal inr the strange 
current was found. On the other hand, Dong and co-workers u'0 used a sum- 
mation of the current insertion over lattice times which requires a further 
identification and fitting of a linear time dependence. They claim to see a 
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Figure 1. Comparison of the electric matrix element calculated with the differential method 
for samples of 100 configurations (open symbols) and 1050 configurations (filled symbols), 
(a) Average of times 10 to 12. (b) Average of times 15 to 17. 



definite signal for the strangeness form factors. 

We present preliminary results from a lattice QCD simulation of the 
strange-quark current loop using a Monte-Carlo sample of gauge field con- 
figurations about ten times larger than in previous studies. Both differential 
and cumulative time analyses are carried out on the simulation data. Com- 
parisons of results obtainpihirom a 100 configuration subsample (the same 
size used in previous workQBi3) with results from the full data set show quite 
clearly how large fluctuations, that could be interpreted as a signal in low- 
statistics data, disappear with improved statistics. Using our complete data 
set consistency between different analysis methods is achieved. By comparing 
different analysis methods and results using different sized gauge field sam- 
ples, criteria for assessing the robustness of a strange quark current signal are 
presented. Using these criteria, no compelling signal for the strange quark 
current is observed. 



2 Lattice Calculations 

The standard methods of the path integral formulation of lattice QCD in 
Euclidean space are used. The Wilson action is used for both quark and gluon 
fields. We need two and three-point functions which describe, respectively, the 
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Figure 2. Electric matrix element using the cumulative method for a sample of 100 config- 
urations. The plots correspond to momentum transfer (a) (1,0,0), (b) (1,1,0), (c) (1,1,1), 
(d) (2,0,0). 



propagation of nucleon states as a function of Euclidean lattice time and the 
strange quark current in the presence of a nucleon. The two-point function 
correlates the excitation of a nucleon state with momentum if at 
some initial time (called 0) and its annihilation at time t. For large Euclidean 
time t this quantity decreases exponentially like e~ Eqt . 

To calculate the three-point function (t,t';Tf ) an insertion of the 
strange-quark vector current is made at time t' . Since there are no strange 
valence quarks in the nucleon this insertion amounts to a correlation of a 
strange-quark current loop with the nucleon two-point function. The strange- 
quark xurrent loop is calculated using a so-called noisy estimator with 
noisetffl. The noise and perturbative subtraction methods employed here are 
the same as those of WilcoxQ except that 60 noises are used instead of 30. 

The three-point function also has an exponential time dependence but is 
more complicated due to the presence of two times. It can be shown that the 
exponential time factors can be cancelled by taking the ratio 
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Figure 3. Electric matrix element using the cumulative method for a sample of 1050 con- 
figurations. The plots correspond to momentum transfer (a) (1,0,0), (b) (1,1,0), (c) (1,1,1), 
(d) (2,0,0). 



m , ^_ g( 3 >(m^)g (2) (*';Q) m 

D- G( 2 )(t;0)G( 2 >(tf;-?) ' 1 ' 

It is these ratios R that are analyzed to get the final results. 

Note that in writing the two and three-point functions, labels associated 
with the Dirac indices of the nuclcon fields and the Lorentz index of the 
current have been suppressed. Also in writing Eq. (|l|) the necessary spin 
sums and projections are not shown, the expression is given schematically to 
show the time dependp&ce of the various factors. The detailed calculations 
follow previous workBtlQ 

The ratio R is usually summed in order to try to improve the signal. The 
differential method uses a difference of R(t, t'; if) on neighbouring time slices. 
It can be shownQ that the quantity M(t, if) given by 

t+i 

M(t, t) = J2 «*) - R tt - ?)) ( 2 ) 

is the current matrix element of interest. 
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Figure 4. Comparison of the magnetic matrix element calculated with the differential 
method for samples of 100 configurations (open symbols) and 1050 configurations (filled 
symbols) using re = 0.152 for the valence quark, (a) Average of times 10 to 12. (b) Average 
of times 15 to 17. 
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Figure 5. Comparison of the magnetic matrix element calculated with the differential 
method for samples of 100 configurations (open symbols) and 1050 configurations (filled 
symbols) using re = 0.154 for the valence quark, (a) Average of times 10 to 12. (b) Average 
of times 15 to 17. 



An alternative is to simply sum R and then fit to a linear time dependence 
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Figure 6. Magnetic matrix element using the cumulative method for a sample of 100 con- 
figurations with k = 0.152 valence quark. The plots correspond to momentum transfer (a) 
(1,0,0), (b) (1,1,0), (c) (1,1,1), (d) (2,0,0). 



to get the matrix element. We use 



t'—t 

S(t,if) = ^ #(*,*';"?), 
t'=i 



constant + tM(t, if). 



(3) 
(4) 



Summijig current insertions up to t' = t follows the suggestion of VichD 
et alS3 and helps to reduce the statistical noise. Dong and co-workersa 
actually use a different upper limit (t' = tuxedit fixed > t). 



3 Results 



Calculations were carried out in quenched approximation at gauge field cou- 
pling of (3 = 6.0. The lattice size was 20 3 x 32. A total of 1050 gauge 
field configurations were generated using a pseudo-heat bath algorithm. Two 
thousand sweeps were used between saved configurations. 

Since quenched lattice QCD does not provide a perfect description of 
hadrons, there is some ambiguity in fixing the parameters (overall scale and 
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Figure 7. Magnetic matrix element using the cumulative method for a sample of 1050 
configurations with k = 0.152 valence quark. The plots correspond to momentum transfer 
(a) (1,0,0), (b) (1,1,0), (c) (1,1,1), (d) (2,0,0). 



quark masses) in the calculations. In this work we use 0.152 as the hop- 
ping parameter for the strange quark. Using a -1 = 2GeV and the </>-meson 
to fix the strange quark mass would suggest a hopping parameter closer J.i 
k s =0.153. On the other hand, using the scale of Dong and co-workersH 
a -1 = \.7AGeV gives k s smaller than 0.152. 

Results for electric matrix element with valence quark k v =0.152 calcu- 
lated using the differential method Eq. (|J) are shown in Fig. (|l|), averaging 
over two different time windows. The 1Q0 configuration sample results are 
consistent with those obtained by WilcoxB. No signal is seen when the gauge 
field sample size is increased to 1050. The summed ratio Eq. (3) is shown 
in Fig. (j^) and ([|) as a function of the nucleon sink time for different mo- 
mentum transfer's. It shows oscillations characteristic of lattice correlation 
function ratiosEl The magnitude of these oscillations decreases slowly as the 
configuration sample size is increased. 

Results for the magnetic matrix element for different valence quark masses 
calculated using the differential method are plotted in Fig. (^) and (^) . As in 
the electric case the results are consistent with zero. 

Finally, the summed ratios for the magnetic current are given in Fig. (^|) 



Presented at EMI 2001, Osaka, Japan 



(a) 



(b) 



Oil 



H h 



(c) 



(d) 



10 15 20 

time 



10 15 20 

time 



Figure 8. Magnetic matrix element using the cumulative method for a sample of 100 con- 
figurations with k = 0.154 valence quark. The plots correspond to momentum transfer (a) 
(1,0,0), (b) (1,1,0), (c) (1,1,1), (d) (2,0,0). 



to Fig. (|). Note that a kinematic factor (see Eq. (3) in Ref.i) of q/(E + M) 
has been removed. The results for 100 configurations, Fig. (||) and (g), should 
be compared to Fig. (1) of Dong et alB and Fig. (2) of Mathur and Dongu 
where calculations with the same statistics are reported. Then, comparing 
with Fig. (0) and (|^), one sees that the. kind of oscillations in the time range 
10 to 15 which Dong and co-workersBQ took to be their signal, have largely 
disappeared in the higher statistics results. Of course, fluctuations still persist 
at larger times but even higher statistics simulations will be necessary to 
establish if they go away or if indeed a strange quark current signal is hiding 
under them. 



4 Conclusions 

To get an estimate of the strange-quark current matrix elements requires the 
extraction of a small signal in the presence of large statistical fluctuations. 
The results presented here suggest a number of useful criteria that should be 
met before one can claim a credible signal: 
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Figure 9. Magnetic matrix element using the cumulative method for a sample of 1050 
configurations with k = 0.154 valence quark. The plots correspond to momentum transfer 
(a) (1,0,0), (b) (1,1,0), (c) (1,1,1), (d) (2,0,0). 



• There should be consistency between different analysis methods. 

• The signal should appear in the same lattice time region and its statistical 
significance should increase as the size of the Monte-Carlo sample of gauge 
fields is increased. 

• The signal should appear in the same time window for different masses. 

Not all of these criteria have been met at our present level of statistics. 
Work is continuing and it is hoped to have final results with an increased 
configuration sample size in the not too distant future. 



Acknowledgments 

This is work is supported in part by the National Science Foundation under 
grant 0070836 and in part by the Natural Sciences and Engineering Research 
Council of Canada. We thank N. Mathur for helpful discussions. 



Presented at EMI 200 1, Osaka, Japan 



11 



References 



1. K. F. Liu, J. Phys. G 27, 511 (2001). 

2. D. H. Beck and R. D. McKeown, arXiv:hcp-ph/0102334 contains an ex- 
tensive review. 

3. R. Hasty et al. [SAMPLE Collaboration], Science 290, 2117 (2000). 

4. K. A. Aniol et al. [HAPPEX Collaboration], Phys. Lett. B 509, 211 
(2001). 

5. W. Wilcox, Nucl. Phys. Proc. Suppl. 94, 319 (2001). 

6. S. J. Dong, K. F. Liu and A. G. Williams, Phys. Rev. D 58, 074504 
(1998). 

7. N. Mathur and S. J. Dong, Nucl. Phys. Proc. Suppl. 94, 311 (2001). 

8. S. J. Dong and K. F. Liu, Nucl. Phys. Proc. Suppl. 26, 353 (1992). 

9. S. J. Dong and K. F. Liu, Phys. Lett. B 328, 130 (1994) 

10. J. Viehoff et al. [SESAM Collaboration], Nucl. Phys. Proc. Suppl. 63, 
269 (1998). 

11. S. Aoki et al. [JLQCD Collaboration], Nucl. Phys. Proc. Suppl. 47, 354 
(1996). 



Presented at EMI 2001, Osaka, Japan 12 



